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ABSTRACT 

Iron,  the  universe’s  most  abundant  refractory  element,  is  highly  depleted  in  both  cireumslcttar  and  interstellar 
environments,  meaning  it  exists  in  solid  form.  The  nature  of  this  solid  is  unknown.  In  this  Letter,  we  provide 
evidence  that  metallic  iron  grains  are  present  around  oxygen-rich  asymptotic  giant  branch  stairs,  where  it  is  obser- 
vationally  manifest  as  a  featureless  mid-infrared  excess.  This  identification  is  made  using  Spitzer  Space  Telescope 
observations  of  evolved  globular  cluster  stars,  where  iron  dust  production  appears  ubiquitous  and  in  some  eases 
can  be  modeled  as  the  only  observed  dust  product.  In  this  context,  FeO  is  examined  as  the  likely  carrier  for  the 
20  pm  feature  observed  in  some  of  these  stars.  Metallic  iron  appears  to  be  an  important  part  of  the  dust  condensation 
sequence  at  low  metallicity,  and  subsequently  plays  an  influential  role  in  the  interstellar  medium.  We  explore  the 
stellar  metallieities  and  luminosities  at  which  iron  formation  is  observed,  and  how  the  presence  of  iron  affeets  the 
outflow  and  its  chemistry.  The  conditions  under  which  iron  can  provide  sufficient  opacity  to  drive  a  wind  remain 
unclear. 

Key  words:  circumstellar  matter  -  globular  clusters:  individual  (NGC  362,  NGC  5139,  NGC  5927)  -  infrared: 
stars  -  stars:  AGB  and  post-AGB  -  stars:  mass-loss  -  stars:  winds,  outflows 


I.  INTRODUCTION 

Iron  is  the  sixth  most  abundant  element  in  the  universe,  and 
the  most  abundant  refractory  element.  It  is  observed  to  be  highly 
depleted  in  both  interstellar  and  eireumstcllar  environments 
(Savage  &  Bohlin  1979;  Maas  et  al.  2005;  Delgado  Inglada 
et  at.  2009),  and  must  therefore  predominantly  exist  in  an 
unknown  solid  form.  Iron  may  be  incorporated  in  other  dust 
grains,  primarily  silicates,  though  these  are  usually  iron-poor 
(Gait  &  Sedtmayr  1999;  Kemper  et  al.  2002).  Alternatively,  iron 
may  form  a  metallic  condensate  (Kemper  et  al.  2002;  Verhoelst 
et  al.  2009). 

Dust  grains  form  around  evolved  stars  by  condensation  from 
the  gas  phase,  either  directly,  or  onto  molecular  “seeds”.  Indi¬ 
vidual  dust  species  can  usually  be  identified  by  distinet  infrared 
emission  bands.  In  oxygen-rieh  environments,  amorphous  or 
crystalline  forms  of  silicates,  spinels,  and  corundum  are  ob¬ 
served  to  form,  while  earbon-rieh  environments  give  rise  to 
amorphous  carbon  (amC),  graphite,  and  SiC.  Metallic  iron 
grains,  however,  produce  featureless  infrared  emission  which 
can  be  difficult  to  differentiate  from  other  sources,  particularly 
amorphous  carbon  dust.  Iron  has  hitherto  only  been  inferred  as 
a  likely  component  of  oxygen-rich  dusty  winds  (Kemper  et  al. 
2002;  Verhoelst  et  al.  2009),  but  never  positively  identified. 

Previous  observational  studies  (McDonald  et  al.  2009;  Boyer 
et  at.  2009)  of  giant  stars  in  globular  clusters  have  found 
a  featureless  contribution  in  addition  to  the  flux  emanating 
from  the  star’s  photosphere,  which  has  been  attributed  to  an 
unidentified  eireumstellar  dust  species.  We  herein  show  that 
this  dust  species  is  metallic  iron. 


2.  EVIDENCE  FOR  DUST  EMISSION 
2.1 .  Observations 

Our  sample  consists  of  35  highly  evolved,  (strongly)  pulsat¬ 
ing  stars  in  globular  clusters  (Sloan  et  al.  2010,  Paper  I),  ob¬ 
served  using  the  Spitzer  Space  Telescope  Infrared  Spectrograph 
(IRS).  These  long-period  variable  stars  lie  near  or  above  the 
red  giant  branch  tips  of  19  Galactic  globular  clusters,  which 
range  in  metallicity  from  IFe/H]  =  -1.62  to  -0. 10.  Ob¬ 
served  with  Spitzer's  low-resolution  modules,  the  spectra  cover 
a  =  5.2-38  pm  with  a  resolving  power  of  R  ~~  60-120. 

One  carbon-rich  star  (Lyngd  7  VI)  is  discounted,  leaving  a 
sample  of  34  stars.  Of  these,  silicate  features  at  10  and  20pm 
are  seen  in  24  objeets,  implying  oxygen-rich  chemistry.  More 
pertinently,  the  remaining  10  stars  either  have  very  weak  silicate 
features  or  appear  to  be  “naked”. 

JHKL-band  light  curves  are  also  available  (Paper  I)  for  25  of 
our  34  objects,  which  allows  us  to  construct  the  photometric 
fluxes  of  these  objeets  at  the  pulsation  phase  of  their  IRS 
spectrum. 

Figure  1  shows  the  L-[8]  color  for  these  objects,  where  the 
8  jrm  photometry  was  derived  by  convolving  the  IRS  spectrum 
with  the  Spitzer  IRAC  8/zm  filter.  A  truly  naked  photosphere 
should  have  L  —  [8)  0,  once  interstellar  reddening  has  been 

taken  into  account7  (cf.  Boyer  et  at.  2008).  All  25  stars,  however, 
have  a  positive  color.  This  implies  a  flux  excess  at  8  pm  which 
cannot  be  eaused  by  eommonty  assumed  dust  species  (silicates. 


7  We  here  assume  E(L  —  [81)  =  0.01 7 E(B  -  V).  afler  Whiiiet  (1992)  and 
Flaherty  el  al.  (2007). 
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Figure  1.  Extinct  ion -corrected  L  —  [8]  colors  in  our  sources.  Plus  signs  show 
stars  with  pure  iron  dust,  filled  circles  show  stars  also  hosting  silicate  dust 
(symbol  size  is  proportional  to  strength  of  emission).  Crosses  show  stars 
exhibiting  a  strong  20 nm  feature.  Color  does  not  vary  significantly  with 
E(B  —  V ),  showing  that  this  is  not  merely  a  reddening  effect. 


alumina)  under  normal  conditions  of  temperature  and  grain  size. 
This  excess  flux  at  long  wavelengths  is  confirmed  by  comparison 
to  marcs  photospherie  model  spectra  (Gustafsson  et  al.  1975, 
2008;  McDonald  et  al.  2009)  as  shown  below. 

From  our  sample,  four  example  stars  have  been  selected  for 
detailed  modeling:  NGC  5927  V3,  a  highly  evolved  asymptotic 
giant  branch  (AGB)  star  with  a  strong,  narrow,  10/xm  silicate 
feature;  NGC  6352  V5,  a  less  luminous  giant  star  with  unusual 
dust  features;  and  NGC  362  V2  and  NGC  5139  V42,  “naked” 
RGB-tip  or  AGB  stars  showing  optieal  TiO  bands  (indicating 
they  are  oxygen-rich)  and  infrared  dust  excesses,  but  no  silicate 
features  (Smith  ct  al.  1999;  McDonald  et  al.  2009;  Boyer  et  al. 
2009). 


2.2.  Fitting  the  Observ  ations 

For  a  dust-free  comparison,  we  use  a  stellar  atmosphere  model 
created  using  the  marcs  code  (Gustafsson  et  al.  1975,  2008)  at 
3500  K,  [Fe/H]  =  -1.0,  [a/Fe]  =  +0.3  (see  McDonald  et  al. 
2009  for  details  of  these  models).  Being  very  similar  objeets  near 
the  giant  branch  tip,  we  expect  the  effective  surface  temperature 
below  the  dust-produeing  zone  to  be  3750  ±  ~250  K.  We 
use  the  near-infrared  flux  to  fit  the  photospherie  eontinuum. 
A  low-temperature,  metal-rieh  model  therefore  represents  an 
intentional  bias  toward  a  model  with  a  more-luminous  infrared 
spectrum  compared  to  its  near-infrared  flux.  This  reduces  the 
calculated  excess  to  rule  out  disparity  between  the  modeled  and 
real  photospherie  spectra. 

Figure  2  shows  the  observed  spectra,  alongside  fits  to  the 
stellar  wind  contribution  using  dusty  (Nenkova  et  al.  1999). 
Using  the  marcs  model  atmosphere  as  an  input,  we  fit  a  stellar 
wind  model  using  metallic  iron  grains  (Ordal  et  al.  1988).  For 
each  model,  we  assume  a  radiatively  driven  wind;  and  a  standard 
Mathis-Rumpl-Nordsieek  grain  size  distribution  (Mathis  et  al. 
1977),  given  by  n(a)  =  a~q ,  q  =  3.5  for  a  =  5-250  nm. 

We  correct  for  interstellar  reddening  using  the  absorption 
profiles  from  McClure  (2009).  We  further  assume  that  A\  = 
3.2 E(B  —  V)  and  that  A\ / A k  =  7.75.  E(B  —  V)  values  for 
NGC  362,  5927,  and  6352  are  from  Harris  (1996),  E(B  ~V)  = 
0.10  is  assumed  for  co  Cen,  based  on  Harris  (1996)  and 


McDonald  et  al.  (2009).  The  correction  applied  is  thus  quite 
small  and  negligibly  different  from  that  applied  in  Paper  I. 

NGC  362  V2  and  co  Cen  V42  are  best  fit  with  a  dust  shell 
truncated  to  an  outer  radius  (r^t)  of  ~  1 00  times  the  inner 
radius  (rjnner  =  the  dust  formation  radius).  (Normally,  fits  are 
extended  to  romec/rjnner  ^  1000.)  This  suppresses  the  flux  at 
A  >  30 /zm.  A  similar  effect  can  also  be  achieved  by  using  a 
grain  distribution  with  fewer  large  grains:  both  effects  may  be 
in  play  (we  only  model  a  truncated  shell  here).  The  strength, 
and  indeed  existence  of,  this  suppression  depends  crucially  on 
the  subtraction  of  the  infrared  background  (caused  by  Galactic 
interstellar  medium  and  other  cluster  objeets).  While  every  effort 
has  been  made  to  make  sure  this  has  been  done  eorreetly,  we 
eannot  be  conclusively  sure  that  such  a  suppression  and  thus 
truncation  exist. 

In  NGC  5927  V3  and  NGC  6352  V5,  we  also  observe 
other  dust  species,  which  emit  at  A  >  8/zm.  Identified  species 
inelude  silicates  (magnesium-rieh  olivines  and  pyroxenes)  at 
10  and  20 (jLm  and  amorphous  alumina  (ALOV)  at  1 1-1 5 /zm 
(optieal  constants  from  Draine  &  Lee  1984;  Begemann  ct  al. 
1997).  Related  crystalline  minerals  with  a  smaller  degree  of 
amorphization  are  likely  responsible  for  the  sharp  9.6  and 
13.1  fxm  features  seen  in  NGC  6352  V5  (Sloan  et  al.  2003). 
Note  that  the  fit  to  NGC  5927  V3  produces  a  fairly  good  match 
to  the  10  fxm  feature,  but  at  a  wavelength  which  is  too  red.  This 
can  be  caused  by  a  number  of  factors,  ineluding  grain  size  and 
porosity  (Voshehinnikov  Sc  Henning  2008).  The  20  gm  feature 
has  an  unknown  carrier:  a  promising  candidate  is  magnesio- 
wiistite  (MgrFei-tO;  see,  e.g.,  Poseh  et  al.  2002;  Lebzelter 
et  al.  2006:  Paper  1).  This  provides  a  good  fit  to  the  feature  in 
NGC  6352  V5.  The  comparatively  red  peak  wavelength  seen 
in  NGC  6352  V5  further  suggests  that  x  %  0;  i.e.,  that  it  is 
near-pure  FeO  (Henning  et  al.  1995).  We  return  to  this  feature 
later. 

We  could  not  fit  a  single-temperature  dust  model  to  these  two 
stars,  implying  that  the  dust  species  are  not  in  thermal  contact 
with  each  other.  Iron  and  alumina  were  fit  at  temperatures  of 
^1000  K.  In  NGC  6352  V5,  silicates  and  FeO  were  fit  at 
^450  K;  in  NGC  5927  V3,  the  strong  10 gm  peak  suggests 
mainly  glassy  silicates,  which  we  fit  at  ~400  K. 

Table  1  lists  our  fitted  parameters,  including  dust  temperatures 
at  the  outer  and  inner  edges  of  the  dust  envelope  (7oUlcr,  Tinner) 
and  the  size  of  the  envelope  used  (r(HJler/ /*jnrKr,  r\nncT  2  /?*).  We 
do  not  inelude  mass-loss  rates,  as  explained  at  the  end  of  this 
section. 


2.3.  Alternative  Sources  of  Excess 

Featureless  infrared  excess  ean  be  produced  by  a  variety 
of  circumstellar  material.  This  can  take  the  form  of  free-free 
emission,  emission  from  shells  of  molecular  gas  (a  “molsphere”; 
e.g.,  Tsuji  2000),  extremely  large  silicate  grains,  or  amorphous 
carbon  grains.  Figure  3  shows  modeled  fits  for  these  possible 
emission  sources,  modeled  in  an  identical  way  to  the  iron  fits 
listed  above. 

Ionization  in  a  chromosphere  or  shock  fronts  propagating 
through  the  outer  atmosphere  can  give  rise  to  free-free  emission, 
as  observed  by  radio  observations  toward  Galaetie  Miras  (Reid 
&  Menten  1997).  The  optical  depth  of  free-free  emission  in  a 
region  of  thickness  r  can  be  calculated  as 


Tv 


3.28  x  1<T7 


-1.35 
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GHz 
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Figure  2.  Modeled  dust  contributions*  (solid  lines)  to  Sparser  IRS  spectra  (black  dots)  and  JHKL  photometry  (black  pluses),  which  has  been  corrected  to  the  pulsation 
phases  of  the  IRS  spectra.  The  dashed  line  is  a  marcs  model  denoting  the  photosphenc  contribution  The  bottom  two  panels  include  the  extra  indicated  contributions 
Irom  AI2O3,  silicates,  and  FeO  (dotted  lines). 


Table  I 

Parameters  of  DUSTY  Models  Presented  in  Figure  2. 


Parameter 

NGC  362  V2 

to  Cen  V42 

NGC  5927  V3 

NGC  6352  V5 

Assumed  luminosity  (Lq) 

1826 

1862 

2000 

2000 

Assumed  [Fe/Hl 

-1.16 

-1.62 

-0.37 

-0.70 

Assumed  E(B  -  V) 

0.05 

0.10 

0.45 

0.21 

References 

1,3 

2 

4 

4 

Wind  component 

A 

B 

A 

B 

Composition 

100% 

100% 

100% 

100% 

96%  Fe 

70%  Silicate 

Iron 

Iron 

Iron 

Silicate 

4%  AI2O3 

30%  FeO 

Tmncr  du>t  envelope  edge  (K) 

800 

1000 

1100 

400 

1000 

450 

Toutcr  du.st  envelope  edge  (K) 

194 

173 

Z outer/  dinner 

50 

120 

1000 

1000 

100 

100 

Optical  depth  al  0.55  /on 

0.28 

0.70 

0.55 

045 

1.05 

0.12 

References,  (1)  Boyer  et  al.  2009;  (2)  McDonald  et  al.  2009;  (3)  E{B  -  V)  from  Harris  1996:  (4)  luminosity  estimated, 
metal licity  from  Harris  1996. 


where  the  electron  density,  ne>  is  recoverable  from  the  electron 
temperature,  Tc ,  via  the  Saha  equation.  At  low  frequencies, 
where  xv  1,  the  emission  scales  as  Fv  oc  v~~  for  an 
isothermal,  tsobaric  region;  and  Fy  a  u-06  for  an  isothermal 
wind  expanding  with  constant  velocity.  For  10  Cen  V42,  we 
find  our  spectrum  would  require  zv  =  1  at  3-4  /im,  and  that 
F,  a  n“16  for  rv  1.  This  more-elosely  approximates  an 
isobaric  density  profile,  which  would  be  difficult  to  achieve  in 
an  expanding  wind  Very  similar  conditions  would  have  to  exist 
in  all  34  of  our  stars,  as  the  shape  of  their  spectra  are  very  similar. 

An  isobaric  chromosphere  producing  a  spectrum  with  this 
optical  depth  cannot  exist  except  in  the  following  case:  the  wind 


must  be  mostly  ioni/cd  (>6000  K)  and  either  the  chromosphere 
extends  to  an  unphysical  radius  (>70/?o)  and/or  the  mass-loss 
rate  becomes  unphysical  (>10-55Afo  yr_l)  for  a  thermally 
expanding  wind.  Departing  from  an  isobaric  state  makes  con¬ 
ditions  even  more  prohibitive.  While  chromospheres  may  be 
present  in  these  stars,  this  infrared  excess  is  not  seen  tn  chromo- 
spherically  active  stars  that  arc  less  evolved  (McDonald  &  van 
Loon  2007;  Boyer  et  al.  2009).  In  V42,  circumstcllar  ioniza¬ 
tion  appears  dominated  by  pulsation  shocks.  Reid  &  Mcnten 
(1997)  show  that,  for  Galactic  Miras,  the  temperature  and 
density  increases  provided  by  shocks  in  the  outer  atmosphere 
are  insufficient  to  produce  an  ionized  layer  that  would  be 
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a)  Cen  V42:  molsphere,  -  1021  cm'2, 1400  K.  3  R. 
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to  Cen  V42: 50  iim  silicate  grains  at  1500  K 
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Figure  3.  Fils  to  lhe  Spitzer  IRS  specimm  of  co  Cen  V42  using  aliemathe  emission  sources.  In  each  panel,  lhe  following  is  shown — dois:  tRS  speclrum;  plus  sign: 
phase-corretied  /.-band  phoiometry;  dolled  line*  fil  for  indicated  grain  composition  or  emission  mechanism 


optically  thick  in  the  infrared.  We  thus  rule  out  free-free  emis¬ 
sion  as  a  dominant  contributor  to  our  speetra. 

The  molsphere  creates  a  dense  “forest”  of  molecular  lines. 
Water  lines  are  the  most  important,  and  cause  emission  through¬ 
out  the  spectrum,  but  mainly  in  the  6-8  jim  and  >  20 /zm 
regions.  This  emission  is  seen  in  some  of  our  example  stars 
(Paper  I).  However,  the  water  bands  in  our  speetra  are  very 
weak  compared  to  stars  with  naked  molspheres  (e.g.,  p  Peg — 
Tsuji  etal.  1997).  Gas-phase  model  emission  speetra  (Figure  3), 
ereated  from  the  line  list  of  Partridge  &  Sehwenke  (1997)  using 
SpeetraFaetory8  (Cami  et  al.  2010),  fail  to  reproduce  the  flatness 
of  the  speetrum  for  any  reasonable  combination  of  temperature 
and  column  density.  Figure  3  shows  the  highest  eolumn  den¬ 
sity  we  consider  feasible,  and  even  this  contains  ~10-7  A/q  of 
water  or  ~0.01%  of  the  star’s  oxygen.  Furthermore,  most  mol¬ 
sphere  host  stars  still  require  an  extra  opacity  source  (such  as 
iron  grains)  providing  flux  at  6-8  gm  (Verhoelst  et  al.  2009). 
Infrared  molecular  emission  should  therefore  not  make  a  signif¬ 
icant  contribution  to  our  speetra. 

Of  the  commonly  considered  dust  speeies  only  metallic  iron 
and  amC  produee  a  featureless  infrared  eontinuum  that  matehes 
the  observed  “naked”  speetra.  Other  dust  speeies  (e.g.,  graphite) 
produee  either  strong  speetral  features  or  inflections  which  do 
not  fit  observed  speetra  Spectral  features  can  be  damped  by 
unusual  conditions,  notably  exceptionally  large  grains  (Hofner 
2008):  these  eannot  simultaneously  reproduce  both  the  exeess  at 
k  <  8  jim  and  k  >  20 /im  at  the  observed  levels  (Figure  3)  and 
are  not  expeeted  in  metal-poor  stars  where  grain  growth  may  be 
hampered  by  a  laek  of  constituent  elements. 


8  hup://spectrafaciory.nei 


Observationally,  amC  can  also  provide  an  acceptable  fit  to 
our  example  stars  over  most  of  the  speetrum.  The  remaining 
discrepancies  may  be  fixable  by  fine-tuning  the  wind  parameters 
(grain  size,  density  distributions,  temperatures,  background 
subtraction,  etc.).  Detection  of  amC  as  dominant  dust  species 
in  an  oxygen-rieh  environment  would  be  unexpected.  In  giant 
stars,  earbon  and  oxygen  are  bound  into  CO  near  the  stellar 
surfaee,  locking  away  the  least  abundant  element.  The  remaining 
majority  element,  either  earbon  or  oxygen,  dominates  the  dust 
chemistry.  CO  may  be  dissociated  in  atmospherie  shock  waves 
or  by  ehromospherie  UV  irradiation  (Beck  etal.  1992).  Thisean 
lead  to  the  production  of  small  amounts  of  other  earbon-rieh 
molecules:  mainly  CO2  (observed  in  our  stars:  see  Paper  I),  but 
also  HCN  and  CS  (Duari  et  al.  1 999).  Small  quantities  of  earbon- 
rieh  dust  may  then  form  (Hofner  &  Andersen  2007),  though 
oxygen-rieh  dust  production  is  still  expeeted  to  be  dominant.  In 
our  sample,  however,  the  featureless  component  often  appears 
to  be  the  dominant  speeies,  and  no  evidence  for  oxygen-rieh 
dust  is  found. 

Having  ruled  out  the  alternatives,  we  conclude  that  the  source 
of  the  dust  exeess  must  be  iron  grains,  and  that  these  must 
dominate  the  eondensed  material  in  some  stars.  If  we  naively 
assume  identical  populations  of  spherieal  grains,  the  fraction  of 
warm  silicates  in  the  “naked”  stars  being  <  1  %  of  the  iron  eontent 
(Figure  3;  the  same  would  be  true  at  a  <3%  level  for  an  amC 
component).  The  silieate  fraction  may  increase  if  iron  grains  are 
made  smaller  and  oblate,  thus  presenting  a  higher  eross  seetion 
per  unit  mass:  extreme  needle-like  grains  are  required  if  silicates 
are  to  become  the  dominant  dust  speeies  by  mass.  As  we  do  not 
know  the  grain  shape,  it  becomes  difficult  to  estimate  a  mass- 
loss  rate  for  an  iron  wind.  Spherical  grains  produee  unfeasibly 
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low  expansion  velocities  of  order  ~m  s~  1  for  a  pure  iron  wind; 
expansion  velocities  for  NGC  5927  V3  and  NGC  6352  V5  were 
not  calculated  as  the  model  assumes  a  zero  initial  velocity  for 
both  temperature  components,  which  will  not  be  the  ease  if 
they  form  at  different  radii.  The  mass-loss  rates  and  expansion 
velocities  can  also  be  significantly  modified  by  altering,  e.g.,  the 
grain  size  or  density  (single  crystal  versus  conglomerate),  but 
reasonable  changes  still  fail  to  produce  feasible  values,  except 
for  extremely  elongated  grains. 

This  implies  that  iron  alone  may  not  be  able  to  drive  a 
dusty  wind  in  these  stars.  Other  factors  may  also  accelerate 
the  wind.  These  could  include  small  amounts  of  amorphous  car¬ 
bon  (Hofner  &  Andersen  2007);  energy  deposition  either  from 
shocks  created  by  stellar  pulsations  (Bowen  1988),  or  magnetic 
reconnection  (Pijpers  &  Hearn  1989);  increased  opacity  from 
eircumstellar  molecular  gas  (Elitzur  et  al.  1989);  stellar  pul¬ 
sation,  either  by  direct  acceleration  or  dissipation  of  shocks; 
or  line-driving  of  H  l  or  Ca  n  atoms,  or  H2O  molecules  (e.g., 
Elitzur  et  al.  1989;  Bowen  1989).  The  driving  mechanism  is 
thus  likely  to  be  complex,  with  radiation  pressure  on  metallic 
iron  being  only  one  of  several  accelerants  which  drive  matter 
from  the  star.  Obtaining  an  accurate  radial  profile  of  the  w  ind's 
outflow  velocity  will  be  crucial  in  determining  whteh  factors 
dominate  this  acceleration. 

3.  DISCUSSION 

3.1.  The  Fate  of  Iron  and  Iron -oxide  Grains 

Efficient  production  of  metallic  iron  by  oxygen-rich  AGB 
stars  would  explain  several  iron  depletions.  Gas-phase  depletion 
of  iron  extends  to  the  interstellar  medium  (Savage  &  Bohlin 
1979).  Iron  grains  contribute  to  infrared  emission  from  this 
medium  (Chlewicki  &  Laurcijs  1988)  and  iron  is  seen  to  be 
destroyed  in  shocks  (Leboutcillcr  et  al.  2008).  Iron  is  also  de¬ 
pleted  by  >90%  in  the  gas  phase  of  planetary  nebulae  (Delgado 
lnglada  et  al.  2009).  Direct  evidence  for  iron  condensation  in 
AGB  winds  comes  from  depletion  patterns  in  binary  post-AGB 
stars.  In  these  stars,  the  wind  is  captured  in  a  eireumbinary 
disk,  from  where  gas  re-accretes  onto  the  star.  These  stars  show 
large  depletions  of  refractory  elements,  with  iron  being  among 
the  most  depleted  (Maas  ct  al.  2005)  thus  indicate  a  history  of 
iron-rich  dust  production. 

The  20  /xm  feature,  seen  in  NGC  6352  V5,  is  also  observed 
in  NGC  5927  VI  and  Terzan  5  V6.  It  has  a  width  (FWHM)  of 
~3/xm  and  a  line-to-eontinuum  ratio  up  to  50%.  This  feature  has 
previously  been  seen  in  low  mass-loss  stars;  it  correlates  with 
the  presence  of  an  1 1  /im  shoulder,  and  13.1  and  28  //m  peaks 
(Littlc-Marenin  &  Little  1988;  Sloan  et  al.  2003).  The  related 
peaks  are  also  present  in  our  sample,  and  we  confirm  that  the 
20  /xm  feature  disappears  for  more  evolved  (i.e.,  more  luminous) 
stars.  The  20  /xm  and  associated  features  are  barely  detectable 
in  the  naked  stars,  strongest  in  stars  with  weak  silicate  emission 
and  disappear  when  silicates  dominate  (Figure  4).  This  suggests 
they  are  produced  either  during  the  formation  or  destruction  of 
silicate  grains.  The  identification  of  the  feature  with  Mg?Fei_vO 
is  unconfirmed  mainly  because  it  does  not  prediet  the  associated 
features,  but  the  proposed  identifications  for  those  (corundum, 
pyroxenes;  Sloan  et  al.  2003)  do  not  provide  good  fits  at  20  /xm, 
making  it  plausible  that  a  range  of  minerals  form  under  these 
conditions.  We  suggest  that  these  features  may  be  eaused  by  the 
chemical  modification  of  silicate  grains:  silicate  grains  which 
adsorb  iron  into  the  lattice  have  higher  opacity,  and  can  be 
heated  until  the  iron-bearing  component  dissociates  from  the 


Figure  4.  Variation  ofdusl  properties  wilh  bolomeiric  magnitude  and  meiallicily 
Symbols  as  shown  in  Figure  1 


grain  ( Woitke  2006).  This  is  expected  to  occur  at  any  radius  from 
the  star  at  which  dust  formation  is  efficient  (Gail  &  Sedlmayr 
1999).  In  this  model,  metallic  iron  is  produced  first,  but  at 
higher  mctallicity  and  mass-loss  rates,  the  iron  becomes  partly 
incorporated  into  other  grain  types  including  iron-oxide  and 
iron-containing  silicates. 

Studies  of  inclusions  within  meteoric  primitive  solar  system 
material  have  traced  dust  grains  to  specific  origins  in  AGB 
stars,  supemovae  and  novae,  through  their  isotopic  ratios  (Zinncr 
2003;  Messenger  ct  al.  2005).  Recently,  the  first  pre-solar  iron- 
oxide  grain  has  been  identified  (Ross  et  al.  2008).  The  iron  and 
oxygen  isotopic  ratios  in  it  are  indicative  of  an  origin  in  an  AGB 
star  This  supports  the  hypothesis  that  AGB  stars  are  a  source  of 
interstellar  iron  dust,  and  have  contributed  to  the  solar  system's 
iron  content. 

3.2.  The  Dust  Condensation  Sequence  at  low  IFe/H] 

Figure  4  shows  the  variation  of  dust  types  (see  Paper  1)  with 
bolometrie  magnitude  and  metallieity.  Silicate-rich  winds  are 
only  seen  at  high  luminosity  or  high  metallieity,  with  iron  dom¬ 
inating  the  outflow  in  lower-luminosity,  lower-metallieity  stars. 
This  indicates  that  the  dust  condensation  sequence  depends  on 
mctallicity  and  luminosity.  Gail  &  Sedlmayr  (1999)  show  that, 
under  equilibrium  conditions,  iron  is  only  expected  to  condense 
before  silicates  in  regions  with  high  gas  pressure.  If  this  applies 
to  our  stars,  it  implies  high  gas  pressure  at  low  metallieity  and/ 
or  low  luminosity,  as  may  be  expected  if  the  outflow  velocities 
where  iron  condenses  (at  the  outside  edge  of  the  molsphere) 
arc  very  low.  Under  non-equilibrium  conditions,  grain  growth 
chemistry  is  strongly  controlled  by  kinetics  (Taehibana  et  al. 
2009),  which  may  favor  the  formation  of  iron. 

In  Gail  &  Sedlmayr  (1999),  silicates  form  before  iron  at  lower 
densities  and  temperatures.  The  restriction  of  silicates  toward 
higher  metallieity  and  luminosity  (Figure  4)  suggests  that  their 
eircumstellar  envelopes  are  more  extended,  favoring  silicates. 
Silicate  production  is  expected  to  increase  with  both  stellar 
metallieity  and  luminosity  (van  Loon  2000). 

4.  CONCLUSIONS 

Using  mid-infrared  spectra  and  photometry,  we  have  shown 
that  there  is  a  considerable  amount  of  unattributed  infrared 
flux  in  a  large  selection  of  globular  cluster  giants.  Through 
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radiative  transfer  modeling,  we  deduce  that  this  flux  is  most 
likely  due  to  metallic  iron  grains  forming  in  a  truncated  stellar 
wind,  corroborated  by  a  potential  identification  of  FeO  with  the 
20  /xm  emission  feature.  The  production  of  metallic  iron  seems 
to  become  more  dominant  at  lower  metallicity,  suggesting  that 
the  dust  condensation  sequence  is  fundamentally  different  in 
metal-poor  objects.  Large-scale  production  of  iron  grains  and 
iron  oxide  in  AGB  stars  can  explain  iron  depletion  in  the  gas 
and  solid  phases  of  the  post-AGB  stars  and  planetary  nebulae; 
as  well  as  isotopic  ratios  in  FeO  grains  in  meteorites.  While 
iron  increases  opacity  in  oxygen-rich  winds,  it  remains  unclear 
whether  it  can  dominate  the  driving  of  metal-poor  winds.  Future 
observations  to  determine  outflow  velocities  in  these  stars  should 
help  determine  the  driver  of  these  stellar  winds. 

This  Letter  uses  observations  made  using  the  Spitzer  Space 
Telescope ,  operated  by  JPL,  California  Institute  of  Technology 
under  NASA  contract  1407  and  supported  by  NASA  through 
JPL  (contract  number  1257184).  We  thank  Martha  Boyer  for 
her  helpful  comments. 
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